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The crystal structure and lattice phonons of (BEDT-TTF)2l3 superconducting /3-phase (where 
BEDT-TTF is bis-ethylen-dithio-tetrathiafulvalene) are computed and analyzed by the Quasi Har- 
monic Lattice Dynamics (QHLD) method. The empirical atom-atom potential is that successfully 
employed for neutral BEDT-TTF and for non superconducting a-(BEDT-TTF) 2 I 3 . Whereas the 
crystal structure and its temperature and pressure dependence are properly reproduced within a rigid 
molecule approximation, this has to be removed account for the specific heat data. Such a mixing 
between lattice and low-frequency mtra-molecular vibrations also yields good agreement with the 
observed Raman and infrared frequencies. From the eigenvectors of the low-frequency phonons we 
calculate the electron-phonon coupling constants due to the modulation of charge transfer (hopping) 
integrals. The charge transfer integrals are evaluated by the extended Hiickel method applied to all 
nearest-neighbor BEDT-TTF pairs in the ab crystal plane. From the averaged electron-phonon cou- 
pling constants and the QHLD phonon density of states we derive the Eliashberg coupling function 
a(u)F(uj), which compares well with that experimentally obtained from point contact spectroscopy. 
The corresponding dimensionless coupling constant A is found to be ~ 0.4. 
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■ I. INTRODUCTION 

> 

The strength and importance of electron-lattice phonon (e-lph) coupling in the superconductivity mechanism 
on : of organic superconductors has always been rather controversial. Early numerical estimates based on simplified 
models gave very low values for the coupling to acoustic phononsp and much more attention was then devoted 
to electron-molecular vibration e-mv cnupling.El On the other hand, "librons" were also invoked in the pairing 
mechanism of organic superconductors On the experimental side, most of the data have been collected for bis- 
ethylen-dithio-tetrathiafulvalene (BEDT-TTF) salts, which are the most extensive and representative class of ocganic 

■ superconductars.cl In particular, recent Raman experiments on BEDT-TTF salts pointed out that the intensity^ and 
the frequencyu of some low- frequency phonon mode change at the superconducting critical temperature TL. Oddly 
enough, also one mfra-molecular BEDT-TTF mode has been shown to exhibit a frequency shift at T c a Carbon 
isotopic substitution on the central double bond BEDT-TTF was claimed to have dramatic effects on the T c of one 
superconducting BEDT-TTF salt, I but subsequent extensive isotopic substitution studies on other superconducting 
BEDT-TTF salts strongly suggested that the lattice phonons are likely involved in the-auperconducting mechanism.E2l 
Attempts to take into account both e-mv and e-lph coupling have been put forwardJU but the role and the relative 
importance of the two types of coupling in the pairing mechanism is far from being settled. i_. _ . 

Whereas extensive studies have been devoted to the characterization of mira-molecular phonons of BEDT-TTFtJaEfl 
and to the estimate of the relevant e-mv coupling strengthJLd very little is known about the lattice phonon structure in 
BEDT-TTF salts or in other organic superconductors. Obtaining a sound characterization of BEDT-TTF salts lattice 
phonons is not easy, since in general the unit cell contains several molecular units, and the phonon modes obviously 
differ for different crystalline-stpictures. We have tackled the problem by adopting the "Quasi Harmonic Lattice 
Dynamics" (QHLD) method,Ej~lI2l by which we are able to analyze both the crystal and the lattice phonon structure 
in terms of empirical atom-atom potentials, in principle transferable among crystals containing the same atoms. We 
have first obtained Gr-S and H atom-atom potential parameters reproducing crystal structure and lattice phonons of 
neutral BEDT-TTF Jl3 Then we have cpnsidered the I 3 salts, which have only one additional atom to parametrize, 
and present several crystalline phases.u After the successful application of the potential to non-superconducting a- 
(BEDT-TTF)2l3 crystal,E3 we present in this paper the results relevant to the extensively studied superconducting 
/3-phases. 
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/9-(BEDT-TTF)2l3 has been the first ambient pressure BEDT-TTF based superconductor to be discovered^ and its 
unit cell contains only one formula unit.Ei The BEDT-TTF radicals are arranged in stacks, and the stacks form sheets 
parallel to the ab crystal plane. The centrosymmetric linear LJ anions separate the sheets, forming an insulating 
layer. Several variants of the /3-(BEDT-TTF)2l3 phase have been reported, making difficult a full and detailed 
characterization. The electrochemically prepared /3-(BEDT-TTF)2i3 exhibits ambient pressure superconductivity at 
T c = 1.3 K (/? L -(BEDT-TTF) 2 I 3 ) or at T c =8.1 (/3 H - or /3*-(BEDT-TTF) 2 I 3 ) depending on the pressure-temperature 
history of the sample. Such XL increase has been attributed to a pressure induced ordering process of the ethylene 
groups of BEDT-TTF cation.E3 In addition, thermal treatment or laser irradiation ofpthe a-phase yields an irreversible 
transformation to a superconducting phase (T c = 8.0 K), named at-(BEDT-TTF)2l3,cl which was claimed to be similar 
to the /3-phase. Oh the other hand, /3-(BEDT-TTF)2l3 can also be prepared by direct chemical oxidation (f3co~ 
(BEDT-TTF) 2 I 3 ),Eil with T c between 7.1 and-^7.8 K. Recent X-ray data confirm that thermally treated a-(BEDT- 
TTF) 2 I 3 is identical to /3 C o-(BEDT-TTF) 2 I 3 P but it is still not clear whether /3 co -(BEDT-TTF) 2 l3 is the same as 
/3y-(BEDT-TTF)2l3: the possibility of non-stpichiometric phases has also been put forward as an alternative to the 
ordering process in causing a T c of about 8 Kj23 

The paper is organized as follows. We first discuss in some detail the methods we have adopted to calculate the 
structure, the phonon dynamics and the e-lph coupling strength of (BEDT-TTF^Lj salts. The results relevant to the 
/3-(BEDT-TTF)2l3 phase are then presented and compared with available experimental data. Finally, the possible 
role of electron-phonon coupling in the pairing mechanism of organic superconductors is briefly discussed. 

II. METHODS 

A. Quasi Harmonic Lattice Dynamics 

The crystal structure at thermodynamic equilibrium of (BEDT-TTF)2l3 salts is computed using Quasi Harmonic 
Lattice Dynamics (QHLD). In QHLDE£rta the Gibbs free energy G(p, T) of the crystal is approximated with the free 
energy of the harmonic phonons calculated at the average lattice structure {Ti = 1): 

G(p, T) = 4> intcr + P V + + k BT ln 

Here, Winter is the total potential energy of the crystal, pV is the pressure- volume term, X)qi w qj/2 is the zero-point 
energy, and the last term is the entropic contribution. The sums are extended to all phonon modes of wavevector q 
and frequency Uqi. Given an initial lattice structure, one computes $i n tor and its second derivatives with respect to the 
displacements of the molecular coordinates. The second derivatives form the dynamical matrix, which is numerically 
diagonalized to obtain the phonon frequencies w q i and the corresponding eigenvectors. The structure as a function 
of p and T is then determined self-consistently by minimizing G(p, T) with respect to lattice parameters, molecular 
positions and orientations. 

In the case of (BEDT-TTF)2l3 salts, and in particular of the (3- phase, the choice of the initial lattice structure 
is somewhat problematic- due to the conformational disorder of the BEDT-TTF molecules. In fact, the X-rays 
structural investigationsEZI indicate that /3-(BEDT-TTF)2l3 at 120 K is disordered with two alternative sites for 
the terminaLC atoms, labeled 9a, 10a (staggered form) and 9b, 10b (eclipsed form). On the other hand, ab-initio 
calculations^ for neutral BEDT-TTF indicate that the "boat" geometry (C% symmetry) is more stable than the 
"planar" geometry (D2 symmetry) by 0.65 kcal/mole. The "chair" distortion [C s symmetry) is slightly more stable 
than the planar molecule, but still less stable than the boat one. The BEDT-TTF+ ion is planar, and in (BEDT- 
TTF)2l3 crystals we have a statistical mixture of neutral and ionized molecules. On the basis of the site symmetry 
constraints, we observe that neutral molecule boat and chair geometries correspond to the Leung's configurations 
9a, 10a and 9b, 10b, respectively^!] Thus the conformational disorder observed in most BEDT-TTF salts is readily 
understood: the energetic cost of deforming the molecules is small with respect to the energy gain among different 
packing arrangements in the crystals. To investigate at least partially the effect of conformational disorder on the 
stability of (BEDT-TTF)2l3 phases, we have performed several calculations starting from different initials molecular 
geometries, as detailed in Section III. 
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B. Potential Model 



We have adopted a pairwise additive inter-molecular potential of the form $j n ter = ^Eranbmfc/''™' + 
A mn exp(— B mn r mn ) — C mn /r„„] where the. sum is extended to all distances r mn between pairs m,n of atoms in 
different molecules. The Ewald's methodEj is used to accelerate the convergence of the Coulombic interactions 
QmQn/ r mn- The atomic charges q m are the PDQ (PS-GVB) results of a recent ab-initio Hartree-Fock calculations, Ell 
and are introduced to model both the neutral and ionized forms of the BEDT-TTF molecule. The parameters A mn , 
B mn and C mn involving C, H and S atoms are taken from our previous calculation of neutral BEDT-TTF.II3 Since 
in the chosen modefia C-H parameters are computed from C-C and H-H parameters via "mixing rules" , the same 
procedure is adopted here for all the interactions between different types of atoms._The iodine parameters have been 
derived from 9,10-diiodoanthracene, and successfully tested on a- (BEDT-TTF) 2 1 3 .E3 The complete atom-atom model 
is given in Table I. 



C. Specific Heat 

The constant volume specific heat as a function of T is computed directly from its statistical mechanics expression 
for a system of phonons: 



As usual in these cases, eq. (2) is evaluated by sampling a large number of q- vectors in the first Brillouin Zone (BZ). 

In our first attempts to compute Cy, we sampled over regular grids in the BZ. We have found that for T < 5 
K the statistical noise was still noticeable even after summing over several thousands of q-vectors; the results were 
dependent on the sample size. At large T, on the contrary, statistical convergence was quite fast. This pathology 
can be attributed to the fact that, due to the exponential factor in eq. (2), only the phonons with w q i < ksT give 
a non- negligible contribution to Cy (T) . For very low T, only the acoustic branches of the phonons with q close to 
zero have sufficiently small frequencies. With a regular grid, only a few of these vectors are sampled, and most of the 
computer time is wasted over regions of the BZ that are already well sampled. 

To obtain accurate statistics at a reasonable cost, we have used a Monte Carlo (random) integration scheme, biased 
to yield a larger sampling probability close to q = 0. For computational simplicity, we have chosen a three-dimensional 
Lorentzian probability distribution, £(q) oc (1 + a | q | ) , where a is a width parameter. The bias is compensated 
by using the reciprocal of the sampling probability as the sample weight. With this scheme most of the computer 
effort is spent in the region q«0, where a denser sampling really matters. By summing over about 2000 q-vectors, 
we have been able to reach a satisfactory statistical convergence, in the whole range between 0.1 and 20 K. At high 
T, the results coincide with those obtained by integrating over a grid. At low T, Cy goes as T 3 , as it should when 
the acoustic modes are properly sampled, and does not fluctuate with the sample size. 
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D. Coupling with low- frequency intra— molecular degrees of freedom 

In most calculations for molecular crystals all mira-molecular degrees of freedom are neglected and the molecules 
are maintained as rigid units. This rigid molecule approximation (RMA) is reasonable for small compact molecules, 
like benzene, where all normal modes have frequencies much higher than those of the lattice phonons. 

Since for both and BEDT-TTF several investigationsEjO'Lj suggest that there are low frequency mtra-molecular 
modes, the validity of RMA for (BEDT-TTF)2l3 appears questionable. Therefore, we have decided to relax the RMA 
and ta-investigate the effects of the mira-molecular degrees of freedom. For this purpose we adopt an exciton-likc 
model.Eia To start with, it is convenient to use a set of molecular coordinates Qi describing translations, rotations and 
internal vibrations of the molecular units in the crystal. To each BEDT-TTF molecule of N = 26 atoms we associate 
the following 3N coordinates: 3 mass-weighted cartesian displacements of the center of mass, 3 inertia-weighted 
rotations about the principal axes of inertia, and 3N — 6 = 72 internal vibrations (the normal modes of the isolated 
BEDT-TTF molecule). The I^~ ion, which is linear, has 3 translations, 2 rotations and 4 internal vibrations. In order 
to compute the phonon frequencies, we need all derivatives d 2 Q/dQ r idQ S j of the total potential <J> with respect to 
all pairs of molecular coordinates Q r i and Q S j . Here r and s label molecules in the crystal, while i and j distinguish 
molecular coordinates. 
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The potential $ is made of intra- and inter-molecular parts, $; n tra and Winter- In the exciton model, the diagonal 
derivatives of <&i n tra potential are taken to coincide with those of an isolated molecule: d 2 $i ntra /dQl i = w^. Here u> r i 
is the frequency of the i-th normal mode of the r-th molecule. All off-diagonal derivatives are zero, which means no 
coupling among different normal modes, and no coupling between normal modes and rigid roto-translations. These 
assumptions are correct for the intea-molecular potential at the harmonic level (by definition). 

The coupling between the molecular coordinates is given by $; n ter- For /3- (BEDT-TTF) 2 13, Winter is described 
by atom-atom and charge-charge interactions, which are both functions only of the interatomic distance. Since the 
distance depends on the cartesian coordinates of the atoms, X ra , the derivatives of $i n tcr can be directly computed 
in terms of the coordinates X ra , and then converted to molecular coordinates Q r f. 



Here a and b label the cartesian coordinates of the atoms in molecules r and s, respectively, and the matrix dX pa /dQ P i 
describes the cartesian displacements which correspond to each molecular coordinate Q P i- The displacements cor- 
responding to rigid translations and rotations of the molecules can be derived by simple geometric arguments. The 
displacements associated to the mira-molecular degrees of freedom are the cartesian eigenvectors of the normal modes 
of the isolated molecule. The atomic displacements, together with the inter-molecular potential model, determine 
the coupling between mire-molecular and lattice modes. We remark that the intera-molecular degrees of freedom are 
taken into account only as far as their effects on the vibrational contribution to the free energy are concerned. No 
attempt to decrease the potential energy by deforming the molecules is done. 



In molecular crystals, m£ra-molecular vibrations are assumed to couple with electrons through modulation of on- 
site energies {e-mv coupling). Lattice phonons are instead expected to modulate mainly the inter-molecular charge 
transfer (CT) integral, t, the corresponding linear e-lph coupling constants being defined as: 



where tKL is the CT integral between neighboring pairs K,L of BEDT-TTF molecules, and where Q qj - is the dimcn- 
sionless normal coordinate for the j— th phonon with wavevector q. By relaxing the RMA, as explained above, the 
distinction between low-frequency intea-molecular modes and lattice modes is at least partially lost. On the other 
hand, e-mv coupling by the low-frequen cy . m olecular modes is expected to be fairly small, as suggested by the calcula- 
tions available for isolated BEDT-TTF Eao Therefore, we have assumed that the calculated low-frequency phonons 
of /3-(BEDT-TTF)2l3, occurring between and about 200 cm^ 1 , are coupled to the CT electrons only through the t 
modulation. 

To evaluate the g(KL; q, j)'s, we have followed a real space approach. Adopting the extended Hiickel method, for 
each pair K, L of BEDT-TTF molecules within the /3-(BEDT-TTF)2l3 crystal we have calculated txh as the variation 
of the HOMO energy in going from the monomer to the dimer. Such aa_approach is known to give t values in nice 
agreement with those calculated by extended basis set ab-initio methodsS txL is calculated for the dimer equilibrium 
geometry within the crystal, as well as for geometries displaced along the QHLD eigenvectors. The various g(KL; q, j) 
are then obtained by numerical differentiation. We have considered only the modulation of the four largest t's, all 
along the ab crystal plane. 

In the case of e-mv coupling the overall electron-phonon coupling strength is generally expressed by the small 
polaron binding energy, E™ v — J2 i 9i/ u u where both gj, the i-th e-mv coupling coristant, and Ui, the corresponding 
reference frequency, are quite naturally taken as independent of the wavevector q.li3 Also in the calculation the e-lph 
coupling we have assumed the optical lattice phonons as dispersionless, and have performed the calculations for the 
q = eigenvectors only. Within this approximation, symmetry arguments show that only the totally symmetric (A g ) 
phonons can be coupled with electrons. Thus, the overall e-lph coupling strength for the j'-th lattice optical phonon, 
can again be expressed by the small polaron binding energy relevant to the j-ih phonon: ej = ^2xl(9kl j/^j)- The 
total coupling strength is then given by E l f p = Y^j e j ■ 

For the three acoustic branches we must of course consider the q dependence of the g's, the coupling constants 
being zero for q = 0. We have then calculated the coupling strength (e" c ) at some representative BZ edges in the 




ab 



(3) 



E. e-lph coupling constants and the Eliashberg function 



g(KL;q,j) = {dtKL/dQ^) 



(4) 
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a*b* reciprocal plane. For each branch, we have averaged the found e" c , and assumed a linear dependence on |q|. The 
latter assumption is correct only in the small |q| limit.El 

The most important single parameter characterizing the strength of elect pan-phonon coupling in the superconduc- 
tivity mechanism is the dimensionless electrqik-phonon coupling constant A£-3'E3 This parameter is in turn related to 
the Eliashberg coupling function 



A = 2 



a 2 (w)F(uj) 



(5) 



where F(uj) is the phonon density of states per unit cell, and a 2 (u>) is an effective coupling function for phonons of 
energy lo. The e-lph Eliashberg coupling function can be evaluated from the QHIJ3 phonon density of states and from 
the electron-phonon matrix element g(k, k'; j) expressed in the reciprocal space:E-3 



a 2 (Lu)F(cu) = N(E F ) ]T(|g(k,k'; j)\ 2 S(lo - Uqj )) FS (6) 
j 

where q = k' — k, k and k' denoting the electronic wavevectors, and N(Ef) is the density of states per spin per unit 
cell at the Fermi level. In eq.(6), ( )fs indicates the average over the Fermi surface. 

We have calculated the g's in real space, as detailed above. In order to introduce the dependence on the electronic 
wavevector k, as required in eq. (6), we have to describe the electronic structure of the p££-phase metal. To get a 
simple yet realistic model we make resort of the rectangular tight-binding dimer model,E3 where the BEDT-TTF 
dimers inside the actual unit cell are taken as a supermolecule. Actually, as in K-phase, in the /3-phase structure 
BEDT-TTF dimers are clearly recognized (in the present formalism, they correspond to the tAB CT integral) . In this 
model there is only one half-filled conduction band in the first BZ, whose dispersion relation as a function of the txL 
CT integrals is easily obtained: 



e(k) = tAB + tAH cos(k x ) +tAE cos(ky) +tAc cos(k x + ky) (7) 

1 potential is 
idence in rec 

is given byO 



The chemical potential is obtained numerically from the half-filling condition. Within our tight-binding approxima- 
tion, the dependence in reciprocal space of the coupling constants associated to the mier-dimer (inter-cell) hoppings 



g(k, k'; j) = 2i g(KL; q, j) [sin(k + q)R - sin kR] (8) 

where R represents the nearest-neighbor lattice vectors (a, b, a + b), and g(KL;q,j) are the three corresponding 
real space inter-cell CT integrals. The Fermi surface average of eq. (6) can now easily performed numerically for the 
inter-dimer contribution. The coupling constants associated with the modulation of t.Tpr. it intra-dimer CT integrals 
are treated as intramolecular coupling constants, and as such are independent of kJldEZI We finally remark that the 
e-mv Eliashberg coupling function is simply given by [a 2 (uj)F(uj)] e - mv — (iV(i?f4/-A0 2* ffl'K^ ~ w i)> N being the 
number of molecules per unit cell and gt being the usual e-mv coupling constant)!^ so that A e _ ml , = N(Ef)E™ v . 



III. RESULTS 



A. Crystallographic structures 



The unit cell^rAXBEDT-TTF^Ia contains one Ig ion at the (0 0) inversion site and two BEDT-TTF molecules 
at generic sites.EBESEH At 4.5 K the two BEDT-TTF molecules have a boat geometry (with the terminal C atoms in 
9a, 10a positions) and are interconverted by the inversion.Ej At 100 and 120 K the lattice is disordered^] and inversion 
symmetry is satisfied only statistically, with a mixture of boat and chair molecules. 

As explained in section II, at first we have made calculations with rigid molecules, and then we have relaxed the 
RMA with the addition of a subset of mtra-molecular modes. The crystal structure is only marginally affected, by 
RMA and in Table II we report the comparison between RMA-calculated and experimental crystal structuretJcJ at 
several temperatures and pressures. Fig. 1 reports a more extensive and direct comparison between calculated and 
experimental crystal axis lengths against T and p. The calculations have been performed by minimizing the free energy 
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G with the molecules kept rigid at their experimental, ordered geometry at 4.5 K.E2I To investigate the effect of small 
changes in molecular geometry, the structures of /?- (BEDT-TTF) 2I3 have been recomputed with the 120 K geometryEZI 
and ordered molecules (staggered or boat form). The effect of the change in molecular geometry is negligible. At all 
temperatures, /?-(BEDT-TTF)2l3 appears to be thermodynamically more stable than a-(BEDT-TTF)2l3,E9 giving 
account for the irreversible interconversions of a-(BEDT-TTF)2l3 into /3-like phases. 

The effect of molecular deformations has been investigated, within the RMA approximation, by testing several 
model geometries in /?-(BEDT-TTF)2l3 as wsjlLas in a-(BEDT-TTF)2l3 phases. The potential energy has been 
minimized with the experimental geometriesLjCj and with chair-a, chair-/?, boat-a|-and boat-/? model geometries. 
The chair-a geometry is the average of the two chair molecules in a- (BEDI-TTF^FsJ^ while chair-/? is the molecule 
of /3-(BEDT-TTF)2l3 with all the terminal carbons in 9b, 10b positions.E3 The boat-a and /? geometries are those 
observed in the corresponding phases.E3!23 For both a and (3 phases, the potential energy minimum is found with 
experimental geometry of that phase, and the system becomes less stable if any other geometry is used. It should 
be noticed that for /3-(BEDT-TTF)2l3 the boat-/? geometry coincides with the experimental geometry at 4.5 K, and 
thus yields the lowest energy. The difference between a-(BEDT-TTF) 2 I 3 and /?-(BEDT-TTF) 2 I 3 essentially vanishes 
if the chair-a and boat-a geometries are used in the /3-phase, while the chair-/? and boat-/? geometries drastically 
destabilize the a-phase. This behavior clearly indicates that molecular deformations play a crucial role in stabilizing 
the various (BEDT-TTF) 2 I 3 phases. 



B. Specific heat 



We next turn our attention to the phonon structure. In the /?-phase we have only one formula unit in the triclinic 
unit cell, and within RMA we expect 8A g and 6A U q = lattice phonons active in Raman and in IR, respectively. 
The number of phonons experimentally observed in the 10-150 cm -1 spectral region is in any case smaller than the 
above prediction, so vibrational spectra do not offer a very stringent test of the calculations. On the other hand, 
there is another observable, the specific heat, which depends on the frequency distribution. As shown in Fig. 2, at 
20 K the Cy calculated within RMA (dotted line) is about 50% smaller than the experimental C p (dots, from Ref. 
|39| ). The difference between Cy and C p is usually small for solids, since their thermal expansion is small. Therefore, 
we attribute most of the discrepancy between Cy and C p to the infra-molecular modes, which are neglected in the 
RMA calculation. i - in 

Ab-initio calculationsHLj indeed indicate the presence of several low- frequency BEDT-TTF infra-molecular (in- 
ternal) normal modes. Since calculations refer to a free molecule, a direct comparison with experimental data in the 
solid state is not feasible. However, they constitute a very convenient starting point for relaxing RMA in QHLD 
calculations, as explained in Section II. We have included the lowest. nine BEDT-TTF internal modes which fall in 
the same spectral region as the lattice modes (below ~ 220 cm _1 ),t3 and therefore are likely coupled. In addition, 
the symmetric and antisymmetric stretchings, and the two bendings of Ig~, expected at 114, 145, 52 and 52 cm -1 , 
respectively,!!!] have been included in the QHLD calculations. The cartesian displacements of BEDT-TTF were ob- 
tained from the ab-initio calculations while those of I3 were determined by symmetry alone, as often it happens 
for small molecules with high symmetry. 

The Cy computed by relaxing the RMA is also shown in Fig. 2. The agreement with experiment is greatly 
improved with reseect to RMA calculations. We anticipate that the same kind of result has been obtained for k- 
(BEDT-TTF)2I3,e3 and conclude by stating that RMA has to be relaxed for a realistic calculation of the low-frequency 
phonons of BEDT-TTF crystals. 



C. Phonon assignments 

We now go back to the characterization of individual low-frequency phonons. In the RMA classification, below 220 
cm -1 we expect Raman activity for 8 lattice modes, 9 BEDT-TTF infra-molecular modes and one stretching of I^~; 
in IR we expect 6 lattice modes, 9 BEDT-TTF infra-molecular modes ancLtiujce I3 modes. The modes calculated at 
the minimum G structure at 120 K are compared with experimental onesj3'E3cj in Tables III and IV for A u and A g 
modes, respectively. We have chosen the 120 K temperature since in this, way we can compare the normal state phonon 
frequencies and eigenvectors for the minimum G and the experimentalEZI structure. The frequency differences between 
minimum G and experimental structure are quite small. The comparison between calculated and experimental 
vibrational frequencies is satisfactory, although not very significant given the low number of observed frequencies. 
Since we also have all the corresponding eigenvectors, we report an approximate description of the phonons, given 
for both BEDT-TTF and I^T as percentage of the lattice (rigid molecule) and of the infra-molecular contributions. 
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Since in some cases there is a considerable mixing between lattice and molecular modest clear distinction cannot 
be made. Fig. 3 reports the full dispersion curves along the C, V, X, and Y directionsEj and density of states of 
/3-(BEDT-TTF)2l3. In order to make the figure more readable, we have limited the highest frequency to 150 cm -1 . 
Fig. 3 puts in evidence the complex structure of /3- (BEDT-TTF) 2I3 low-frequency phonons. We have a very dense 
grouping of modes in the 50-100 cm -1 region, with several avoided crossings between the dispersion curves, and clear 
mixing between lattice and molecular modes. Only the acoustic phonon branches contribute to the density of states 
below ~ 25 cm -1 , so that the typical uJ 2 dependence is observed. On the other hand, at energies higher than ~ 140 
cm -1 the almost dispersionless mira-molecular modes dominate, and the phonon density of states appears as a sum 
of delta- like peaks (not shown in the figure). 

D. Electron-phonon coupling 

The e-lph coupling constants for the optical phonons of /3-(BEDT-TTF)2i3 are reported in Table IV. As explained 
in Section II, if one assumes that the eigenvectors are independent of q, only the A g phonons can couple to electrons. 
In Table IV for each phonon we report both the individual g(KL; j, q) (Eq. 5) and the small polaron binding energy 
ej. The two most strongly coupled modes are those calculated at 32 and 113 cm -1 . Whereas the former mode has 
been observed in the Raman spectrum, and together with the lower frequency mode (22_cm _1 ) undergoes a drastic 
intensity weakening at T c ,u the latter has not been reported even in the normal stateflo The reason might be due 
to the proximity of the very intense, resonantly enhanced band at 121 cm -1 , due to the symmetric stretch of the 
Ig anion. One band at 107 cm -1 has been observed below 6 K for 488 nm laser excitation.cfl On the other hand, 
a band at 109 cm _1 p-,whose intensity varies with sample and irradiation, has been attributed to the splitting of the 
I3 stretching modejla as a consequence of the commensurate superstructure reported in one X-ray investigation at 
100K.EI Certainly the 100-130 cm -1 spectral region deserves further experimental scrutiny with the latest generation 
of Raman spectrometers. A second observation is that whereas the 113 cm -1 mode involves only the BEDT-TTF 
units, and is mostly a lattice mode, the 32 cm -1 one is a mixing between rigid Ijj" motion and "flexible" BEDT-TTF 
vibrations. This finding suggests a not marginal role of the counter-ions sheets in /3-type BEDT-TTF salts. 

As shown by Table IV, the coupling of individual optical modes with electrons is in general not particularly strong, 
but on the whole the strength of e-lph coupling, as measured by the sum of the e,-, is appreciable, around 45 meV, to 
which we have to add the contribution of the acoustic phonons. For the sake of comparison, we give also the e" c for 
the three acoustic branches, calculated as average over several points at the BZ edges. In order of decreasing phonon 
frequency (see Fig. 3) the e oc are: 2.3, 3.3 and 18.2 meV, respectively. The coupling strength of the lowest acoustic 
branch at the zone edge is comparable to that of the most strongly coupled optical phonons. Thus, the overall e-lph 
coupling strength is of the same order of magnitude as that due to e-mv coupling, about 70 meVJia 

We can make a more direct connection with superconducting properties by calculating the Eliashberg function 
and the dimensionless electron-phonon coupling constant A. As seen in eq. (6), the absolute value of the Eliashberg 
function depends on the electronic density of states at the Fermi energy, N(Ep). Experimental estimates of this 
critical parameter are problematic, since the measured quantities already include or the A enhancement factor, or the 
Coulomb enhancement factor, or both. The available theoretical estimates are all based on the extended Hiickcl tight 
binding method. The choice of N(Ep) = 2.1 spin states/eV/unit cell, as obtained by this methodjla is consisten with 
our extended Hiickel estimates of the CT integrals in real space. To our advantage, we can compare the calculated 
a(uj)F(oj) with that derived from normal state current/voltage measurements at a point contact junction.laej This 
kind of experiment is rather difficult to perform on organic crystals like /3-(BEDT-TTF)2l3, since one has to be careful 
about pressure effects at the point contact. The use of a contact between two /3-(BEDT-TTF)2l3 crystals made the 
current- voltage, characteristics rather stable, from which the a(aj)F(uj) function reported in the upper part of Fig. 4 
was obtained.c3 We have changed the scale on the ordinate axis to maintain the same energy unit (cm -1 ) throughout. 
It is clear that the spectral resolution of the experiment is larger than ^10 cm -1 , and probably increases with 
energy. Indeed, no spectral detail is visible beyond 240 cm -1 , where the contribution of e-mv coupled modes should 
be detectable.B Therefore, to make easier the visual comparison with the experimental data, we have smoothed the 
calculated a(uj)F(u>) (Fig. 4, lower part) by a convolution with a Gaussian distribution. We have also assumed that 
the Gaussian distribution width increases linearly with uj (from 0.1 to 20 cm" 1 in the 1-200 cm -1 interval). 

Fig. 4 puts in evidence the very good agreement between experiment and calculation. The absolute scale of 
a(uj)F{u)) turns out to be practically the same, even if both experiment and calculation are affected by considerable 
uncertainties, as explained above. The three main peaks observed in the experiment are well reproduced and are 
identified as due to the most strongly coupled phonon branches, namely, the optical phonons at 113 and 32 cm , 
and the lowest frequency acoustic branch. The calculated peak frequency due to the latter is slightly higher than the 
experimental one (22 vs 10 cm -1 ). This discrepancy might be due to the fact that the experiment refers to the @l- 
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(BEDT-TTF) 2 I 3 phase, whereas our calculation refers to a perfectly ordered phase like the /?*-(BEDT-TTF) 2 I 3 @ We 
also remark that, at variance with traditional superconductors, the Eliashberg function is remarkably different from 
the phonon density of states (Fig. 3). For instance, the peak around 120 cm -1 in F(u>) is due to the dispersionless 
Eg" stretching mode, which is completely decoupled from the electron system, whereas the broad peak in a(u))F(u>) 
is due to the nearby (113 cm -1 "lattice" mode of the BEDT-TTF molecules. Due to the complex phonon structure, 
a(u>) is not nearly constant, but varies rapidly with the frequency. 

The dimensionless coupling constant A obtained by integration of a(u))F(u))/u> up to 240 cm -1 turns out to be 
around 0.4. The contribution to XJxpm e-mv coupled modesEj is instead around 0.1. Thus in the McMillan picture 
the overall A of /3-(BEDT-TTF) 2 I 3 B is - 0.5, which may well account for the observed T c = 8.1 KE 

IV. DISCUSSION AND CONCLUSIONS 

The computational methods we have adopted to analyze the crystal and lattice phonon structure, and the electron- 
phonon coupling strength of /3-(BEDT-TTF)2l 3 are empirical or semiempirical. The form of QHLD atom-atom 
potentials has no rigorous theoretical justification, aad the corresponding parameters are derived from empirical 
fittings. We have adopted ab-initio atomic chargesEj to take into account Coulomb interactions between atoms, 
and ab-initio vibrational eigenvectorstil to introduce the coupling between lattice and molecular modes. Also the 
extended Hiickel method used to characterize the (3- (BEDT-TTF) 2 I 3 electronic structure is semiempirical, albeit with 
an experimental basis far wider than QHLD. In view of the obvious limitations of empirical or semiempirical methods, 
the success achieved in the case of /3-(BEDT-TTF)2l 3 is even beyond our expectations, also considering that none of 
the empirical parameters has been adjusted to fit /3-(BEDT-TTF)2l 3 experimental data. 

Indeed, all the available /3-(BEDT-TTF)2i 3 experimental data have been accounted for. The crystal structure, and 
its variation with temperature and pressure, is correctly reproduced (Fig. 1 and Table II). Useful hints about the 
relative thermodynamic stability of (BEDT-TTF)2l 3 a and /3 phases have been obtained, as well as some indications 
on the effect of BEDT-TTF conformation of the phases stability. The specific heat (Fig. 2) and the few detected 
Raman and infrared bands have been accounted for by including the coupling with low- frequency molecular vibrations. 
Finally, the point contact Eliashberg spectral function has been satisfactorily reproduced (Fig. 4). 

Despite the success, it is wise to keep in mind the QHLD limitations. First and foremost, conformational disorder 
in the crystal structure is not included. This is not a limitation of the QHLD method only, but it is a serious 
one particularly for /3-(BEDT-TTF)2l 3 salts, where disorder plays an important role even in the superconducting 
properties. It is in fact believed that a fully ordered structure is. at the origin of the higher T c (8.1 K) displayed by 
/?*- or (3h- phases with respect to /3 L (BEDT-TTF) 2 I 3 (1.5 K)£3 Furthermore, even if the QHLD method is able to 
follow the T and p dependence of the crystal structure, phase transitions implying subtle structural changes may be 
beyond its present capabilities, even for fully ordered structures. The relative stabilities of the phases can indeed be 
reproduced only at a qualitative level. For what concerns the electron-phonon coupling, one has to keep in mind that 
it depends on phonon eigenvectors, and these are obviously more prone to inaccuracies than the energies. Finally, 
extension to other BEDT-TTF salts with counter-ions different from Ig is not obvious, requiring additional atom-atom 
parameters. 

Once the above necessary words of caution about the method are spelled out, we can underline what in any case have 
learned from the present QHLD calculations. So far, in the lack of any description, no matter how approximate, of the 
phonons modulating the CT integrals, only speculative discussions about their role in the superconductivity could be 
put forward, catching at best only part of the correct picture. One of the most important indications coming out from 
the present paper is the need of relaxing the RMA. So from one side-we cannot try to focus on the isolated molecule 
mira-molecular vibrations presumably modulating the CT integral,E3 and on the other side that the "librations" of 
the rigid moleculesu lack of precise meaning. In other words, there is no simple or intuitive picture of the phonons 
modulating the CT integrals. Our results, and the overall mode mixing, suggest that also the counter-ions vibrations 
may play a perhaps indirect role in the coupling. 

The results of the present paper definitely assess the very important role played by the low-frequency phonons in 
the superconducting properties of BEDT-TTF salts. Both acoustic and optic modes modulating the CT integral are 
involved. The overall dimensionless coupling constant is ~ 0.4, much larger then that due to it e mv coupled phonons 
(~ 0.1). Of course, a mere numerical comparison of the two A's is not particularly significant, since one has to keep in 
mind the very different time scales (frequencies) of the two types of phonons. The phonons appreciably modulating 
the CT integrals fall in the 0-120 cm -1 spectual region (Table IV), whereas those modulating on-site energies have 
frequencies ranging from.j4.00 to 1500 cm -1 £3 Applicability of the Migdal theorem to the latter appears dubious: 
non-adiabatic correctionscJ or alternative mechanisms such as polaron narrowingc£l have been suggested. For these 
reasons we will not get involved into detailed discussions about the relative role of e-lph and e-mv coupling in the 
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superconductivity mechanism. We limit ourselves to state that phonon mediated coupling can well account for the 
observed critical temperature of the ordered A|(BEDT-TTF)2l3 phase, given plausible values of the other fundamental 
parameter, the Coulomb pseudopotential /i.tllca 

The results of the present paper suggest that phonon mediated mechanism is responsible for the superconductivity 
of BEDT-TTF-based salts—The same conclusion was reached on the basis of the solution of the BCS gap equation for 
k phase BEDT-TTF salts.E3 On the other hand, evidences are also accumulating towards non-conventional coupling 
mechanisms in organic superconductors, such as spin- fluctuation mediated superconductivity.cj Similar apparently 
contrasting experimental evidences are also found for cuprates,E3 pointing to a superconductivity mechanism where 
both electron-phonon coupling and antiferromagnetic spin correlations are taken into account. 
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TABLE I. Parameters for the atom-atom potential Vjj(r) = Aij exp(— Bijr) — dj/r e . V is in 
kcal/mol, r in A, A, B, C in consistent units. Heteroatom parameters are given by Aij — -J ' AnAjj, 



Bij 


= (Sii + B jj )/2 and Cij 


— \J CaCjj. 




i 




Ba 




H 


2868 


3.74 


40.2 


C 


71460 


3.60 


449.3 


S 


329600 


3.31 


5392.0 


I 


642376 


3.09 


16482.6 



TABLE II. Structural data of /3-(BEDT-TTF) 2 i3 as a function of T (K) and p (GPa): experimental (Ref. M and computed 
unit cell axis a, b, c (A), angles a, /3, 7 (degrees) and volume V (A 3 ). The lattice is triclinic, space group PI (Cj ), with Z = 1. 



T 


P 




a 


b 


c 


a 


P 


7 


V 


G(p,T) 


4.5 





expt. 


6.519 


8.920 


15.052 


95.32 


96.09 


110.44 


807.6 








calc. 


6.571 


9.147 


15.073 


93.94 


95.07 


111.85 


832.6 


-229.034 


20 





expt. 


6.543 


8.968 


15.114 


95.34 


96.05 


110.30 


819.1 








calc. 


6.571 


9.147 


15.074 


93.94 


95.07 


111.85 


832.6 


-229.056 


120 





expt. 


6.561 


9.013 


15.173 


95.07 


95.93 


110.28 


829.2 








calc. 


6.577 


9.159 


15.095 


93.85 


95.07 


111.82 


835.9 


-231.677 


298 





expt. 


6.615 


9.100 


15.286 


94.38 


95.59 


109.78 


855.9 








calc. 


6.591 


9.189 


15.145 


93.64 


95.08 


111.72 


844.2 


-242.678 


4.5 


1.5 


expt. 


6.449 


8.986 


15.034 


94.79 


96.57 


111.29 


799.1 








calc. 


6.556 


9.117 


15.048 


94.05 


95.14 


111.85 


826.1 


-211.930 


6.1 


4.6 


expt. 


6.433 


8.947 


14.927 


95.15 


96.77 


111.40 


786.1 








calc. 


6.532 


9.069 


15.007 


94.20 


95.27 


111.86 


816.0 


-175.319 



TABLE III. Low energy A u phonons of /3-(BEDT-TTF) + i3 . 



Expt. a 


Calc. 






Approximate description 




cm -1 


cm" 1 




I 3 " (%) 




BEDT-TTF+ 


(%) 






lattice 




internal 


lattice 


internal 




216 










98 




173 










90 




151 








26 


68 


133 


148 






86 




10 


130 


136 






8 




88 


124 


123 
114 








48 
8 


48 

86 


95 


112 








24 


72 


91 


90 








56 


40 




84 


8 




13 


34 


44 




74 






56 


26 


12 


71 


69 








22 


78 




67 


8 




40 


18 


36 




62 


22 




44 


10 


20 




57 






23 


44 


28 


48 


43 


55 




11 


26 






33 


44 






22 


30 




17 


41 






22 


36 



a From Ref. 
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TABLE IV. Low enerj 


;y A g phonons 


and coupling constants of /3-(BEDT-TTF) + I^. 






Expt. a 


Calc. 


Approximate description 




Coupling Constants (meV) b 




ej, meV 


cm -1 


cm -1 


T - (<VA 

h K /o ) 


bedt-ttf; 


r (%) 


r,( A R- n\ 

9\ At >>3) 


^ AC- A\ 


9V Atj i J) 










lattice internal 


lattice internal 












214 


207 






on 
yu 




— i 


-i 

— i 






179 


172 






no 




— i 


— i 








165 






64 


6 




z 


-l 




149 


144 




1 R 


oil 


— 1 


o 


-1 

1 


-2 






131 




1 o 




— 






1 




121 


120 


on 
89 






— 3 


1 

1 


-1 

1 


-1 






113 




~.Q 
•JO 


•i o 


1 Q 


I 


o 
O 


-1 


16 




111 






9U 


— 4 


1 

1 


1 

— 1 


3 




91 


87 


i n 
1U 


50 


36 


o 
I 




3 


4 






81 




48 


48 


— 1 


1 

— 1 


— 1 


-5 






73 




8 


88 


2 




2 


1 






64 




78 


22 


-4 


6 


-3 


3 


9 




60 


33 


18 


50 


1 


-3 


-4 


5 


7 


53 


51 


25 


50 


22 


4 




-6 


3 


10 




49 


20 


66 


14 


4 


-3 


-1 


2 




39 


44 


33 


60 




-1 


-2 


-2 


2 




32 


32 


56 


16 


24 


6 


-2 


-5 


1 


17 


27 


29 


10 


76 


10 


-4 


4 


2 


2 


11 



a From Ref. |,[l§ 

b The modulated hopping integrals between dimers are labeled according to Ref. |%| the calculated equilibrium values are: 
t A B = 0.22 eV, t A c = 0.08 eV, t A E = 0.10 eV, t AH = 0.06 eV. 
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FIG. 1. Calculated and experimental of /3-(BEDT-TTF)2l3 crystallographies axis lengths as functions of temperature and 
pressure. The experimental points are taken from Ref. ^ (diamonds) and from Refs. p4|p^ (asterisks). 
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FIG. 4. Upper panel: the Eliashberg function as measured from point-contact tunneling experiments (adapted from Ref. 
3). Lower panel: the calculated contribution to a 2 {uj)F{uj) from low frequency e-lph coupled phonons. 
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